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Abstract. - We present a comprehensive theory of the magnetization relaxation in a M1112- 
acetate crystal based on thermally assisted spin tunneling induced by quartic anisotropy and 
weak transverse magnetic fields. The overall relaxation rate as function of the magnetic field 
is calculated and shown to agree well with data including all resonance peaks. The Lorentzian 
shape of the resonances is also in good agreement with recent data. A generalized master 
equation including resonances is derived and solved exactly. It is shown that many transition 
paths with comparable weight exist that contribute to the relaxation process. Previously 
unknown spin-phonon coupling constants are calculated explicitly. 



The molecular magnet Mni2-acetate — a spin 10 system with large easy-axis anisotropy — 
has attracted much recent interest |l], |[ || since several experiments on the magnetization 
relaxation have revealed pronounced peaks [Q, ^, |(| in response to a varying magnetic field 
H z applied along the easy axis of the crystal. These peaks occur at field values where spin 
states become pairwise degenerate. Following earlier suggestions ||, this phenomenon has 
been interpreted as a manifestation of resonant tunneling of the magnetization. However, 
in a critical comparison between model calculations |], [ll], [n], |l2], [l4| and experimental 
data|^, |j| [ll| Friedman et a/.[^5) point out that a consistent explanation of the relaxation 
is still missing. One of the main challenges for theory is to explain the overall shape of the 
relaxation curve as well as the Lorentzian shape of the measured resonance peaks |l5|]. 

In this work we shall present a model calculation of the magnetization relaxation which is 
based on thermally assisted tunneling. We find for the first time reasonably good agreement 
both with the overall relaxation rate (including the resonances) measured by Thomas et al. j^j 
(see fig.l) and with the Lorentzian shape of the first resonance peaks (see figs. 2, 3) measured 
by Friedman et a/. |]l5| for four different temperatures. 

The model introduced below contains five independent parameters: three anisotropy con- 
stants A ^> B 3> £?4, a misalignment angle 6 (angle between field direction and easy axis), 
and sound velocity c. Previously unknown spin-phonon coupling constants are determined 
for the first time; we find that, quite remarkably, they can be expressed in terms of A only. 
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The constants A, B, B± arc known within some experimental uncertainty from independent 
measurements Jl6|, is typically about l°|lf|, while c is not measured yet (to our knowl- 
edge). We achieve optimal agreement between our theory and data as follows. In accordance 
with pj] we set 9 = 1°, while the constants A, B,Bi are fitted to the relaxation data by 
observing, however, the constraints that A, B, B4 are allowed to vary only in the range of their 
experimental uncertainties. The sound velocity c, on the other hand, is a free fit parameter, for 
which we obtain c = (1.45 — 2.0) ■ 10 3 m/s, which is within expected range. Thus, in contrast to 
previous calculations ||, [l^, [ll], |l2] , O , p"8| our theory is in reasonably good agreement not only 
with the relaxation data of refs. |5j, |l5| but with all experimental parameter values known so far. 
In addition, a new prediction is made (sound velocity c) which can be tested experimentally. 

Extending previous workj^, |l^, [ll], |l^, [l3| [l8| we make use of a generalized master equation 
which treats phonon-induced spin transitions between nearest and next-nearest energy levels as 
well as resonant tunneling due to quartic anisotropies and transverse fields on the same footing. 
We solve the derived master equation by exact diagonalization. In addition, we identify for the 
first time the dominant transition paths (see fig. 4) and show that the magnetization reversal 
is not dominated by just one single term but rather by several paths of comparable weight. 
Finally we note that our model is minimal in the sense that it is sufficient to explain the 
existing data. Accordingly, we do not include hyperfine fields Jl^, since the width of 



the resulting GaussiansjllJ 12 turns out to be smaller for T > IK than the width of the 
Lorentzians obtained below (cf. [|l9)). Similarly, we ignore dipolar (HI |l8| interactions since 
they have been ruled out by experiments on diluted samples [pp| . 

Model. To describe the behavior of the Mni2-acetate molecule we follow earlier work[|| 
[l0| , [ll| , [12] , |l3| , 21 and adopt a single-spin Hamiltonian which includes spin-phonon coupling. 
Explicitly, H = H a + Hz + Ti- sp + Ht, where 

H & = -AS 2 Z - BSt (1) 

describes the magnetic anisotropy with easy axis (A ^> B > 0) along the z-direction. S 
is the spin operator with s = 10, and A/k B = 0.52 - 0.56 Kg 0, and B/k B = (1.1 - 
1.3) • 10~ 3 K|il| [l7| are the anisotropy constants (ks is the Boltzmann constant). The 
Zeeman term through which the external magnetic field H z couples to the spin S is given 
by Ti-z = —g^BH z S Z i while the tunneling between ^-states is described by 

H T = ~B 4 (SX + St) -gn B H x S x , (2) 

where H x — |H| sin^ (<C H z ) is the transverse field, with 8 being the misalignment angle. Here, 
B 4 /k B =_(4.3 - 14.4) • 10~ 5 K[p|, and g = 1.9||. Finally, the most general spin-phonon 
coupling |23| is described by 

W sp = 9i{txx-t V y)®{Sl- Sy) + )^g2e X y®{S x ,Sy\ (3) 

+ \gz{txz ® {S x , S z } + e yz «) {S y , S z }) + ^g A {w xz ® {S x , S z } + oj yz ® {S y , S z }) , 

where g$ are the spin-phonon coupling constants, and e a ^ (u> a fi) is the (anti-)symmetric part 
of the strain tensor defined through the displacement u as (du a /dxp ± dup/dx a )/2. To 
determine ^ occurring in (^) we follow ref.p3 (full details will be given elsewhere p5[ ) . From 
the displacement resulting from rotation only, u = 6(p x x (in leading order), we obtain 
the infinitesimal rotation angle 5(f) = x u = {u> yz ,u zx ,uj xy ). Rotating then the spin 
vector S we find (to leading order in ui a p) that the easy axis term, —AS Z , is transformed into 
A(ui xz {S x , S z }+ujy Z {Sy, S z }). Comparison with the last term in (Eh then yields 34 = 2A. Next, 
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expanding the rotation matrices and u up to second-order, we find[p5[, —S<fi x = e yy +e zz —e xx , 
and cyclic. After rotating the rhs of —AS 2 = — A(S 2 — S x — S 2 ) we obtain a term of the form 
A (e xx — e yy ) (S x — Sy) . Comparing this with (||) we see that g\ = A, and thus g\ = g^/2 = A. 
Finally, we note that the terms cx 31,2 produce second-order transitions, while the ones oc 33,4 
produce first-order transitions. Thus, it is very plausible, in accordance with to adopt 
the approximations, \g2\ ~ gi = A, and | (73 1 ss g^ = 2 A (the sign is irrelevant for the transition 
rates calculated below). 

We denote by m), — s < to < s, the eigenstate of the unperturbed Hamiltonian 7i a + TLz 
with eigenvalue e m = —Am 2 — Bm A — gpBH z m. If the external magnetic field H z is increased 
one gets doubly degenerate spin states whenever a level m coincides with a level to' on the 
opposite side of the well (separated by the barrier w A). The resonance condition for double 
degeneracy, i.e. e m — e m ', leads to the resonance field 

gp B H™ n ' = -{m + m')[A + B (to 2 + m' 2 )] . (4) 

As usual, we refer to to + to,' = even (odd) as even (odd) resonances. 

Master equation. We describe the relaxation of the magnetization in terms of a master 
equation for the reduced density matrix p(t) which includes off-diagonal terms due to reso- 
nances. Using the notation p mm ' = (to \p\ to'), p m = (to \p\ to), we start from the generalized 
master equation 

Pm'm — ^m'm ^ ^ ^VmnPn 7m' 'mPm' ra ) (5) 

where in the white-noise approximation 7 m ' m = J mm > = (W m + W m >)/2, with W m = 
Ylnt^rr,) W nm |25| . The phonon induced transition rates W mn are calculated below. Eq. 
(|) holds for correlation times r c much smaller than the spin relaxation time r - being 
satisfied here 28|. Next we incorporate resonant tunneling into (|5|). Let |to) and |to') be 



two unperturbed eigenstates on the lhs and rhs of the barrier, resp., which are degenerate 
when SH Z = H™ m — H z vanishes, but get detuned otherwise. In the presence of tunneling, 
induced by Ht, the two states form (anti-) symmetric levels split by E mm / (for SH Z = 0). 
E m m' is obtained by standard means [|29| . The 2-state Hamiltonian including a bias field then 



reads 



Ht = Cm \m) (m\ H — — — \m) (m'\ + (to <-> m'), (6) 



with Cm = gnB&H z m. Ht provides a valid description as long as the level splitting A = 
V (Cm - Cm') 2 + E mm' remains smaller than |e m — £ TO ±i|, and \s m i — e m /±x|- Next we consider 
the quantum dynamics of the associated density matrix adopting standard arguments p7|. 
From the von Neumann equation p — i [p, Ht] /h, we get p m = iE m m' {Pmm' — Pm'm) /2h, 

jimm' +"/mm'JPmm'^ ^ (p m — Pm' ), (7) 

where Cmm' = Cm — Cm', and likewise for to «-> to'. Here we have allowed for damping due 
to phonons (i.e. finite life-time of the levels to, to') according to (||). Inserting the stationary 
solution of (^) into the p m -equation, we get p m = T™ (p m t — p m ), where 

r m' p2 W m + Wm' 

rn mm' , /TT r TT7 - x 9 / , V / 

1 h 2 (Wm + Wm') /4 



~>m III' 



is the transition rate from m to to' (induced by tunneling) in the presence of phonon-damping. 
We note that V™ has a Lorentzian shape with respect to the external magnetic field 8H Z (the 
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independence of W m around the resonances turns out to be much weaker); it is this T™ that 
will determine the peak shape of the magnetization resonances (see below and figs. 1-3) [|30|. 

To obtain the complete master equation we include also the transitions involving non- 
resonant levels and get 

Pm = -W m p m + ^ W mnPn + T™ (p m > ~ Pm), (9) 

and the same with to <-> to' (see rcf. [^5| for details). A similar equation has been used in 
rcf.jio]], whose r™ has been obtained through a different and rather lengthy derivation filf . 
For levels k ^ m, to', eq. @ reduces to p k = -WkPk + J2n w knPn- 

As in ref.p^| we evaluate the spin-phonon rates via Fermi golden rule (with thermal 
averaging over phonons) and fmd[p5| 

TJ7 _ A2s ±l (£m±l-£m) 3 , 1fV , 

Wm±1 - m ~ 12^SV(^i-^)-l' (10) 

W - !7^ 2 S±2 (gm±2 - £m) 3 , . 

m±2 ' m ~ 192tt P cW efi(* m ±*-'n>) - l ' 

where s±i = (s=Fto)(s±to+1)(2to±1) 2 , and s±2 — (s=Fto)(s±to+1)(s=Fto— 1)(s±to+2). The 
mass density p is given by 1.83 ■ 10 3 kg/m 3 |32]]. Here, c is the sound velocity of the Mn-crystal, 
which is the only free parameter in our theory. It has not been considered before and we are 
not aware of any measurements of c. 

Relaxation time. We solve now the master equation exactly to find the largest (i.e. dom- 
inant) relaxation time. To this end it is convenient to rewrite (^) as p(t) — Wp(t), where 
p n are the components of p. The rate matrix W has 21 eigenvalues Wi found by exact 
diagonalization. The eigenvalues Wi turn out to be non-degenerate with the smallest one 
being far separated from the remaining ones. The relaxation time is then obtained from 
r = maxi {— 1/Rc w;}. The result is plotted in fig.l, where the overall relaxation rate r 
is shown as function of H z and at T = 1.9K. We note that in our model the even reso- 
nances are induced by the quartic i?4-anisotropy, whereas the odd resonances are induced 
by product-combinations of B^S±- and iJ^^-terms. In accordance with the experimental 
uncertainty|l5| we keep —1° < 9 < 1° for all plots. Then the relaxation time r at an even 
resonance peak is about 3 times smaller than the one at a subsequent odd resonance peak. 
The relevant tunnel splitting energies of the even and odd resonances are about the same: 
£4,-4 w E2-3 f» -E3.-5 w -4 f» hs ■ 20mK. (^2.-2/^5 ~ 0.5K.) For comparison we also 
include in fig.l the data reported by Thomas et al.^, [54|. We have optimized the fit (as 
explained in the introduction) in such a way that the fits of the model parameters, given by 
A/k B = 0.54K, B/k B = 1.1 • 1CT 3 K, and B 4 /k B = 4.3 • 10~ 5 K, = 1°, are roughly within 



the reported experimental uncertainties) 16 17]] (see above), while the fit of the sound velocity 
yields c = 1.45 • 10 3 m/s. Almost identical plots are obtained for 0.5° < 9 < 3° g ||. 

In figs. 2, 3 we plot the peaks of the first resonance at H z = (induced only by the B4- 
anisotropy) for four different temperatures, T = 2.5, 2.6, 2.7, 2.8K. The peaks (like all others) 
are of single Lorentzian shape as a result of the 2-state transition rate given in (||) (see also 
[p0|). For comparison we plot in figs. 2, 3 the data reported by Friedman et aZ. |l^] for the same 
temperatures (no error bars, however, are given). The optimal fit values are: A/k B = 0.56K, 
B/k B = 1.3 • 1CT 3 K, B 4 /k B = 10.2 • 10" 5 K, for all -1° < 6 < 1°, and c = 2.0 • 10 3 m/s. Note 
that these values are the same for all four temperatures. These fitting parameters turn out to 
be somewhat larger than in fig.l, which could be caused by differences e.g. in volume-to-surface 
ratio and/or in shape anisotropy, etc. Indeed, the sample of ref.[fL5| consists of many small 
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crystallites in contrast to the single crystal used in ref.(5). In an y case, the differences are 
small, and the sound velocity c seems to be within the expected order of magnitude. Clearly, 
it would be highly desirable to check this prediction by an independent measurement of c. To 
sum up, the agreement between theory and experiment is reasonably good, and our model and 
its evaluation contains the essential physics responsible for the magnetization relaxation. 

Relaxation paths. It is instructive to determine the dominant transition paths via which the 
spin can relax into its ground state. For this we derive an approximate analytic expression for 
the relaxation time denoted by r* (to distinguish it from the exact r obtained above). First, 
we solve the master equation for one particular transition path n which does not intersect with 
other paths. For H z > we find (details will be given elsewhere |p5[) : 

X " (12) 



" 1 + C0C".-"- •) Ira ' 

where j m = W mm i , or T™ , depending on the particular path n characterized by the sum over 
the levels m (see fig. 4). Eq. ( |l2| ) holds for arbitrary initial (£i) and final (sf) energies; for 
H z = (i.e. Ei = Sf) eq. (Ill) reduces to the result found by Villain et aZ.jl3]. If there is more 
than one dominant path (which is typically the case) we have to account for intersections at 
vertices. For this we associate with each path a probability current J n = p n , and interprete 



eq. (12) as a serial circuit with the summands playing the role of "resistances". This allows 
us then to set up flow diagrams for J n (see fig. 4), which obey the analog of Kirchhoff's rules: 
(Kl) J n = 0: The sum over all incoming and outgoing currents vanishes at a vertex 
(current conservation). (K2) J2 n Jni~ n = AN: The sum over all voltage drops (J n T n ) is equal 
to the source-drain voltage AN = p s — p- s for any closed path (probability conservation). 
The total probability current is given by J = AN. Fig. 4 shows the complete (a) and its 
serially reduced (b) flow diagram for < H z < A/gps- From (Kl) we get J — J 1 + J 2 , 
•h + J5 = J&, Ji = J3 + -hi J3 + -h = hi J4 = -h + h, and J 7 + J 8 = J, while from (K2) 
we get AN = Jit ± + J 3 r 3 + J 7 r 7 , J 3 r 3 = J 4 t 4 + J 5 r 5 + J 6 r 6 , J 2 r 2 = Jm + J4T4 + J 5 r 5 , 
JsTs = J5T5 + J§tq + J7T7. From these equations we obtain the relaxation time defined by 
r* = AiV/«7p5|. Finally, when plotted as function of H z there is no visible difference between 
the exact r and this approximate r*, which confirms that the diagram in fig. 4 contains the 
physically relevant relaxation paths close to the first resonance. Similar results are obtained 
for the other resonances p5||. 

We are grateful to A. Chiolero and T. Pohjola for useful discussions. This work has been 
supported by the Swiss NSF. 

Additional Remark. After submission of this work, we were notified by J.R. Friedman that 
our prediction of the sound velocity c is in very good agreement with recent specific heat 
measurements p5[ (see [E5j for details). 



REFERENCES 

[1] C. Paulsen et al, J. Magn. Magn. Mater., 140-144 (1995) 379. 

[2] R. Sessoli et al, Nature (London), 365 (1993) 141. 

[3] M.A. Novak et al, J. Magn. Magn. Mater., 146 (1995) 211. 

[4] J.R. Friedman et al, Phys. Rev. Lett, 76 (1996) 3830. 

[5] L. Thomas et al, Letters to Nature, 383 (1996) 145. 

[6] J.M. Hernandez et al, Europhys. Lett, 35 (1996) 301. 

[7] B. Barbara et al, J. Magn. Magn. Mater., 140-144 (1995) 1825. 

[8] M.A. Novak, R. Sessoli, in Quantum tunneling of Magnetization, eds. L. Gunther, B. Barbara 
(Kluwer, Dordrecht, 1995). 



6 EUROPHYSICS LETTERS 

[9] J. Villain et al, Europhys. Lett, 27 (1994) 159. 

[10] A. Fort et al, Phys. Rev. Lett., 80 (1998) 612. 

[11] F. Luis et al, Phys. Rev. B, 57 (1998) 505. 

[12] F. Hartmann-Boutron et al, Mod. Phys. B, 10 (1996) 2577. 

[13] D.A. Garanin, E.M. Chudnovsky, Phys. Rev. B, 56 (1997) 11102. 

[14] L. Gunther, Europhys. Lett., 39 (1997) 1. 

[15] J.R. Friedman et al, Phys. Rev. B, 58 (1998) R14729. 

[16] A.L. Barra et al, Phys. Rev. B, 56 (1997) 8192. 

[17] Y. Zhong et al, condmat/9809133. 

[18] F. Fernandez et al, J. Appl Phys., 83 (1998) 6940. 

[19] W. Wernsdorfer et al, preprint. 

[20] R. Sessoli, Mol. Cryst. Liq. Cryst., 274 (1995) 145. 

[21] J.M. Hernandez et al, Phys. Rev. B, 55 (1997) 5858. 

[22] R. Sessoli et al, J. Am. Chem. Soc, 115 (1993) 1804. 

[23] E. Callan, H. Callan, Phys. Rev., 139A (1965) 455. 

[24] V. Dohm, P. Fulde, Z. Phys. B, 21 (1975) 369. 

[25] M.N. Leuenberger, D. Loss, to be published. 

[26] A. Abragam, A. Bleany, Electron Paramagnetic Resonance of Transition Ions, Clarendon 
Press, Oxford, 1970. 

[27] K. Blum, Density Matrix Theory and Applications, 2nd edition, (Plenum Press, 1996), ch. 8. 

[28] An estimate for thermal phonons yields r c ~ h/ksT ~ 10~ n s, at T = IK, whereas r ~ 1 s. 

[29] D.A. Garanin, J. Phys. A, 24 (1991) L61. 

[30] Note that in the t(H z ) plot these Lorentzians are truncated by the spin-phonon transition rates 
W m and W m i in such a way that the effective linewidth is much larger than (W m + W m >)/2^zE^. 

[31] J. Villain et al, J. Phys. I France, 7 (1997) 1583. 

[32] T. Lis, Acta Crystallogr. Sec. B, 36 (1980) 2042. 

[33] The present theory holds for \H X \ < 1000 G (which is well satisfied here); otherwise the shift of 

the levels \m) due to H X S X must be taken into account |25|| . 

[34] In [^| the correct temperature is T = 1.9K, not 2. IK (B. Barbara, W. Wernsdorfer, priv. comm.). 

[35] A.M. Gomes et al, Phys. Rev. B, 57 (1998) 5021. 



M. N. Loucnbcrgcr et al. SPIN RELAXATION IN MNi 2 -ACETATE 



7 




Fig. 1. - Full line: semilogarithmic plot of calculated relaxation time r as function of magnetic 
field H z at T = 1.9K. The optimal fit values are: A/k B = 0.54K, B/k B = 1.1 ■ 1(T 3 K, and 
B 4 /k B = 4.3 • 1(T 5 K, 9 = 1°, and c = 1.45 • 10 3 m/s. Dots and error bars: data§. 
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Fig. 2. - Full line: plot of calculated relaxation rate T = 1/r as function of H z for the first 
resonance peak at T — 2.6K. The Lorentzian shape is due to T™ in (||). The optimal fit 
values are: A/k B = 0.56K, B/k B = 1.3 • 1CT 3 K, B 4 /k B = 10.2 • 1(T 5 K, -1° < 6 < 1°, and 
c = 2.0 • 10 3 m/s. Dots: data@. 
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Fig. 3. - Full lines: semilogarithmic plots of calculated relaxation rate T = 1/r as function 
of H z for the first resonance peak at (a) T = 2.5K, (b) T = 2.6K, (c) T = 2.7K, and (d) 
T = 2.8K. The optimal fit values are the same as in fig. 2. Dots: datafl5|j. 




Fig. 4. - (a) Spin relaxation paths (from s = 10 to s = —10) for < gfJ,sH z < A. Dashed lines: 
tunneling transitions due to B<±- and iJ x -terms. (b) Associated serially reduced diagram with 
relaxation times r ra given in (H). For \H Z \ < 0.05T only the path n — ► T3 — > T7 is relevant. 



